vignettes/My Examples/Initialize_Dowling_Balanced_Constrained_NaturalStimuli_Example.R

In michaelkriegsman/DiDiSTATIS: Run DiDiSTATIS and related methods

#Initialize_Balanced_Dowling_ArousalValence.R

SortConstrained_36 <- read.csv('C:/Users/Michael A. Kriegsman/Google Drive/Projects/DiDiSTATIS_Composers/WithMoreSubj_ForDissertation/natural36_constrained_forR.csv', row.names=1)
Sort <- SortConstrained_36[-37,]
MusExp <- t(SortConstrained_36[37,]) #subjects' musical experience
colnames(MusExp) <- "MusExp"


#Balance the groups
#Randomly select 7 participants to remove from the Mid Musical Experience group
#Leaves 10 participants in each group.

# drop_these_MidMusExp <- sample(which(DESIGN$tables$MusExp$mat[,2]==1), 7)
#The results were c(3, 22,  9, 21, 31, 24, 11)

# 4.25.2018, I had previously removed subjects 33 and 34, as mentioned below, but these are in the High Experience group, so I need to keep them, and instead I revert my balancing to the previous setup...
# #Then, noticed that bc033 and bc034 each created 15 categories, and their PFS are unstable and unusual and huge.
# #So, take those out, and bring in 2 other subjects
# drop_these_MidMusExp <- c(9, 21, 31, 24, 11, 33, 34)
drop_these_MidMusExp <- c(3, 22,  9, 21, 31, 24, 11)
Sort_balanced <- Sort[,-drop_these_MidMusExp]


#Create the DESIGN matrix on the rows (arousal_valence)
AV <- c(3, 6, 8, 11, 12, 13, 18, 26, 29, 30, 34, 36)
av <- c(2, 7, 9, 25, 28, 32)
Av <- c(10, 23, 24, 27, 31, 33)
aV <- c(1, 4, 5, 14, 15, 16, 17, 19, 20, 21, 22, 35)

Arousal_Valence_DESIGN_mat <- matrix(0, nrow(Sort_balanced), 4, dimnames = list(rownames(Sort_balanced), c('AV', 'av', 'Av', 'aV')))
Arousal_Valence_DESIGN_mat[AV,1] <- 1
Arousal_Valence_DESIGN_mat[av,2] <- 1
Arousal_Valence_DESIGN_mat[Av,3] <- 1
Arousal_Valence_DESIGN_mat[aV,4] <- 1
# rowSums(Arousal_Valence_DESIGN_mat)




#Create the DESIGN matrix on the tables (still Mus.Exp)
MusExp_balanced <- MusExp[-drop_these_MidMusExp]
MusExp_labels <- c('Low', 'Mid', 'High')
#cut the DESIGN variable
MusExp_vec <- cut(x = MusExp_balanced, breaks = c(-1,.5, 5, 50), labels = MusExp_labels) #0-0.9; 1-4.9, 5+
#create the DESIGN matrix, with columns in order
MusExp_mat <- makeNominalData(as.matrix(MusExp_vec))[, paste0(".", MusExp_labels)]
rownames(MusExp_mat) <- paste0("s", 1:nrow(MusExp_mat))
colnames(MusExp_mat) <- MusExp_labels






#Create the DESIGNs!
DESIGN <- Initialize_DESIGNs(DESIGN_rows_mat   = Arousal_Valence_DESIGN_mat,
                             DESIGN_tables_mat = MusExp_mat, colors_D = c('#fcd305','#ed961c','#ce5050'))